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We propose a flexible stochastic framework for modeling the market share dynamics over time 
in a multiple markets setting, where firms interact within and between markets. Firms undergo 
stochastic idiosyncratic shocks, which contract their shares, and compete to consolidate their 
position by acquiring new ones in both the market where they operate and in new markets. 
The model parameters can meaningfully account for phenomena such as barriers to entry and 
exit, fixed and sunk costs, costs of expanding to new sectors with different technologies and 
competitive advantage among firms. The construction is obtained in a Bayesian framework by 
means of a collection of nonparametric hierarchical mixtures, which induce the dependence 
between markets and provide a generalization of the Blackwell-MacQueen Polya urn scheme, 
which in turn is used to generate a partially exchangeable dynamical particle system. A Markov 
Chain Monte Carlo algorithm is provided for simulating trajectories of the system, by means of 
which we perform a simulation study for transitions to different economic regimes. Moreover, it 
is shown that the infinite-dimensional properties of the system, when appropriately transformed 
and rescaled, are those of a collection of interacting Fleming- Viot diffusions. 

Keywords: Bayesian nonparametrics; Gibbs sampler; interacting Fleming- Viot processes; 
interacting Polya urns; market dynamics; particle system; species sampling models 

1. Introduction 

The idea of explaining firm dynamics by means of a stochastic model for the market evolu- 
tion has been present in the literature for a long time. However, only recently, firm-specific 
stochastic elements have been introduced to generate the dynamics. Jovanovic [26] was 
the first to formulate an equilibrium model where stochastic shocks are drawn from a 
distribution with known variance and firm-specific mean, thus determining selection of 
the most efficient. Later Ericson and Pakes [12] provide a stochastic model for industry 
behavior which allows for heterogeneity and idiosyncratic shocks, where firms invest and 
the stochastic outcome determines the firm's success, thus accounting for a selection pro- 
cess which can lead to the firm's exit from the market. Hopenhayn [24] performs steady 
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state analysis of a dynamic stochastic model which allows for entry, exit and hetero- 
geneity. In [38] a stochastic model for market share dynamics, based on simple random 
walks, is introduced. The common feature of this non-exhaustive list is that, despite the 
mentioned models being inter-temporal and stochastic, the analysis and the explicit de- 
scription of the model dynamics arc essentially done at equilibrium, thus projecting the 
whole construction onto a static dimension and accounting for time somehow implicitly. 
Indeed, the researcher usually finds herself before the choice between a dynamic model 
with a representative agent and a steady-state analysis of an equilibrium model with 
heterogeneity. Furthermore, relevant for our discussion are two technical difficulties with 
reference to devising stochastic models for market share dynamics: the interdependence 
of market shares, and the fact that the distribution of the size of shocks to each firm's 
share is likely to depend on that firm's current share. As stated in [38], these together 
imply that an appropriate model might be one in which the distribution of shocks to 
each firm's share is conditioned on the full vector of market shares in the current period. 

The urge to overcome these problems from an aggregate perspective, while retaining the 
micro dynamics, has lead to a recent tendency of borrowing ideas from statistical physics 
for modeling certain problems in economics and finance. A particularly useful example 
of these tools is given by interacting particle systems, which are arbitrary-dimensional 
models describing the dynamic interaction of several variables (or particles) . These allow 
for heterogeneity and idiosyncratic stochastic features but still permit a relatively easy 
investigation of the aggregate system properties. In other words, the macroscopic behav- 
ior of the system is derived from the microscopic random interactions of the economic 
agents, and these techniques allow us to keep track of the whole tree of outcomes in an 
inter-temporal framework. A recent example of such an approach is given in [5], where 
interacting particle systems are used to model the propagation of financial distress in a 
network of firms. Another example is [36], which studies limit theorems for the process 
of empirical measures of an economic model driven by a large system of agents that in- 
teract locally by means of mechanisms similar to what, in population genetics, are called 
mutation and recombination. 

Here we propose a Bayesian nonparametric approach for modeling market share dy- 
namics by constructing a stochastic model with interacting particles which allows us 
to overcome the above mentioned technical difficulties. In particular, a nonparametric 
approach allows us to avoid any unnecessary assumption on the distributional form of 
the involved quantities, while a Bayesian approach naturally incorporates probabilistic 
clustering of objects and features conditional predictive structures, easily admitting the 
representation of agents' interactions based on the current individual status. Thus, with 
respect to the literature on market share dynamics, we model time explicitly, instead 
of analyzing the system at equilibrium, while retaining heterogeneity and conditioning 
on the full vector of market shares. And despite the different scope, with respect to the 
particles approach in [36] , we instead consider many subsystems with interactions among 
each other and thus obtain a vector of dependent continuous-time processes. In con- 
structing the model, the emphasis will be on generality and flexibility, which necessarily 
implies a certain degree of stylization of the dynamics. However, this allows the model 
to be easily adapted to represent diverse applied frameworks, such as, for example, pop- 
ulation genetics, by appropriately specifying the corresponding relevant parameters. As 
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a matter of fact, we will follow the market share motivation throughout the paper, with 
the parallel intent of favoring intuition behind the stochastic mechanisms. A completely 
micro-founded economic application will be provided in a follow-up paper [32]. However, 
besides the construction, the present paper includes an asymptotic distributional result 
which shows weak convergence of the aggregate system to a collection of dependent 
diffusion processes. This is a result of independent mathematical interest, relevant, in 
particular, for the population genetics literature, where our construction can be seen as 
a countable approximation of a system of Fleming- Viot diffusions with mutation, selec- 
tion and migration (see [6]). Appendix A includes some basic material on Fleming- Viot 
processes. 

Finally, it is worth mentioning that our approach is also allied to recent developments 
in the Bayesian nonparametric literature: although structurally different, our model has 
a natural interpretation within this field as belonging to the class of dependent processes, 
an important line of research initiated in the seminal papers of [30, 31]. Among others, 
we mention interesting dependent models developed in [9-11, 35, 39] where one can 
find applications to epidemiology, survival analysis and functional data analysis. See the 
monograph [23] for a recent review of the discipline. Although powerful and flexible, 
Bayesian nonparametric methods have not yet been extensively exploited for economic 
applications. Among the contributions to date, we mention [22, 27, 33] for financial time 
series, [19, 21] for volatility estimation, [28] for option pricing and [3, 8] for discrete 
choice models, [20] for stochastic frontier models. With respect to this literature, the 
proposed construction can be seen as a dynamic partially exchangeable array, so that 
the dependence is meant both with respect to time and in terms of a vector of random 
probability measures. 

To be more specific, we introduce a flexible stochastic model for describing the time 
dynamics of the market concentration in several interacting, self-regulated markets. A po- 
tentially infinite number of companies operate in those markets where they have a positive 
share. Firms can enter and exit a market, and expand or contract their share in competi- 
tion with other firms by means of endogenous stochastic idiosyncratic shocks. The model 
parameters allow for barriers to entry and exit, costs of expansion in new markets (e.g., 
technology conversion costs), sunk costs and different mechanisms of competitive advan- 
tage. The construction is achieved by first defining an appropriate collection of dependent 
nonparametric hierarchical models and deriving a related system of interacting general- 
ized Polya urn schemes. This underlying Bayesian framework is detailed in Section 2. 
The collection of hierarchies induces the dependence between markets and allows us to 
construct, in Section 3, a dynamic system, which is driven by means of Gibbs sampling 
techniques [18] and describes how companies interact among one another within and be- 
tween markets over time. These undergo stochastic idiosyncratic shocks that lower their 
current share and compete to increment it. An appropriate set of parameters regulates 
the mechanisms through which firms acquire and lose shares and determines the compet- 
itive selection in terms of relative strengths as functions of their current position in the 
market and, possibly, the current market configuration as a whole. For example, shocks 
can be set to be random in general but deterministic when a firm crosses upwards some 
fixed threshold, meaning that some antitrust authority has fixed an upper bound on the 
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market percentage which can be controlled by a single firm, which is thus forced away 
from the dominant position. The competitive advantage allows for a great degree of flexi- 
bility, involving a functional form with very weak assumptions. In Section 4 the dynamic 
system is then mapped into a measure- valued process, which pools together the local in- 
formation and describes the evolution of the aggregate markets. The system is then shown 
to converge in distribution, under certain conditions and after appropriate rescaling, to 
a system of dependent diffusion processes, each with values in the space of probability 
measures, known as interacting Fleming-Viot diffusions. In Section 5 two algorithms 
which generate sample paths of the system are presented, corresponding to competitive 
advantage directly or implicitly modeled. A simulation study is then performed to explore 
dynamically different economic scenarios with several choices of the model parameters, 
investigating the effects of changes in the market characteristics on the economic dynam- 
ics. Particular attention is devoted to transitions of economic regimes as dependent on 
specific features of the market, on regulations imposed by the policy maker or on the 
interaction with other markets with different structural properties. Finally, Appendix A 
briefly recalls some background material on Gibbs sampling, Fleming- Viot processes and 
interacting Fleming-Viot processes, while all proofs are deferred to Appendix B. 



2. The underlying framework 

In this section we define a collection of dependent nonparametric hierarchical models, 
which will allow a dynamic representation of the market's interaction. 

Let a be a finite non-null measure on a complete and separable space X endowed with 
its Borcl sigma algebra ^T, and consider the Polya urn for a continuum of colors, which 
represents a fundamental tool in many constructions of Bayesian nonparametric models. 
This is such that X\ ~ a(-) / a(K) , and, for n > 2, 

A„ Ai,...,A„_i ~ — — — , {I) 

a(X) + n — 1 

where 5 y denotes a point mass at y. We will denote the joint law of a sequence 
(Xl, . . .,X n ) from (1) with M%, so that 

" a(X)W a(X) + i-l 1 ' 

In [2] it is shown that this prediction scheme is closely related to the Dirichlet process 
prior, introduced by [16]. A random probability measure P on (X, is said to be a 
Dirichlet process with parameter measure a, henceforth denoted P ~ ^(-|a), if for every 
k > 1 and every measurable partition Bi, . . . ,Bk of X, the vector (P(Bi), . . . , P(Bf.)) has 
Dirichlet distribution with parameters (a(Bi), . . . ,a(Bk)). 

Among the various generalizations of the Polya urn scheme (1) present in the literature, 
a recent extension given in [37] will be particularly useful for our construction. Consider, 
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for every n > 1, the joint distribution 

n 

q n (dx 1 ,...,dx n ) ocp n (dxi,...,dx n ) J\0n{%k), (3) 

fc=i 

where (3 n is a given bounded measurable function on X. A representation for (3) can 
be provided in terms of a Dirichlet process mixture model [29]. In particular, it can be 
easily seen that when p n = .M" in (3), the predictive distribution for Xi, given X(_i) = 

\X\ , . . . j Xi—i , X%-\-\, . . . , %n) ; ^ 

n 

q n<i (dx t |x(_j) ) oc P n (xi)a{dxi ) + ^ j3 n (xi )6 Xk (dxi ) . (4) 

This can be thought of as a weighted version of (1), which is recovered when /3„ = 1 for 
all n > 1. A more general version of (4) can be obtained by making (3 n depend on the 
whole vector and thus allowing for a broad range of interpretations. See the discussion 
following (18) for this and for a more detailed interpretation of (4) in the context of the 
present paper. 

Consider now the following setting. For each n, let (Xi, . . . ,X n ) <G X" be an n-sized 
sample from and let 

n 

a Xu ..., x „(dy) = a(dy) + ^25 Xk (dy). (5) 
fe=i 

Define the double hierarchy 

X 1 ,...,X n \P iA ^ P, P~#(-|a), 

(6) 

Y\ j . . . , Y n \Q<n ~ Qm Qn ^ \ol X \ .. . .,x n ) ■ 

Here (X\, . . . ,X n ) are drawn from a Dirichlet process P ~ 2){-\a) and (Yi, . . . ,Y n ), given 
(Xi, . . . ,X n ), are drawn from a Dirichlet process Q n := Q|(Ai, . . . ,X n ) ~ ^(■|ax 1 ,...,x„)- 
It can be easily seen that the joint law of (Yi, . . . ,Y n ) conditional on (Xi, . . . ,X n ) is 
M.n x x " , with A4™ as in (2). The following result, stated here for ease of reference, 
can be found in [41]. 

Lemma 2.1. Let Ai" be as in (2). Then 

J AC 1 *"(dy 1 ,...,dy n )MZ(dx 1 ,...,dx n )=MZ(dy 1 ,...,dyn)- (7) 

In particular, Lemma 2.1 yields a certain symmetry in (6), so that we could also 
state that the joint law of {X\, . . .,X„) conditional on (Yi, . . ■ ,Y n ) is Mn vl " " Vn ■ Denote 
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x = (x\ , . . . , x n ) and extend (3) to 

n n 

g„(dx) ocp„(dx) Y[P (dy) (xp n 

k=l k=l 

From (4) , when (X\ , X n ) and ( Y\ , . . . , Y n ) come from (6) we have for 1 < i < n 

n 

92n, l (dx i |x ( _ l ),y) oc Pn{xi)a yi ,..., yn (dxj) + J~] /3 n {x k )6 Xk (cLcj) (8) 

and similarly for y^. It is now straightforward to iterate the above argument and allow 
for an arbitrary number of dependent hierarchies. Denote x r = (x\, . . . , x^) and a x r = 
ot-x\,...,x r n i where r,r',r" belong to some finite index set X, whose cardinality is denoted 
#1. Then, for every n > 1, let 

X r |P i ~ d -P ) P~0(-|a), 

X r \P rU ~-p r , P r ~ &{-\a^), 

„ , „■■* , „ , ,, (9) 

xr | pr . iT . i.J^.p r , r ^ ^@(.\a xr ^,), 



where the dimension subscript n has been suppressed in X r , X r , X r , . . . , for notational 
simplicity. Denote now with 

D n = n ■ #Z (10) 

the total number of components in (9). The joint law of the D n items in (9) can be 
written 

M%(dx r )MZ xr (dx r ')M? ^' (dx r ") • • • , (11) 

where, in view of Lemma 2.1, (11) is invariant with respect to the order of r, r',r", 

With a slight abuse of notation, define 

l(- r )={x r ': r 'el,rVr}, (12) 

to be the set of all system components without the vector x r , and 

!(-<) = K': r'6l}\K} (13) 

to be the set of all system components without the item x[. Analogously to (8) in this 
enlarged framework, the predictive law for x\ , conditional on the rest of the system, can 
be written 

g Dnti (d<|I(-<))cx/3„(4)a I( _ r) (d<) + ^^(4)^(d<), (14) 

k^i 
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where the interpretation of cn x( _ r) is clear from (5) and (12). Note that this predictive 
law reduces to (1) when /3 n = 1 and oi x ,_. = a. Expression (14) will be the key for 
the definition of the market dynamics by means of an interacting system of particles. 
A detailed interpretation for qD n ,i will be provided in the following section. See (18) and 
the following discussion. 

To conclude the section, it is worth noting that (9) generates a partially exchangeable 
array, where partial exchangeability is intended in the sense of de Finetti (see, e.g., [4]). 
That is, if r,r',r" identify rows, then the system components are row- wise exchangeable 
but not exchangeable. 

3. Dynamic models for market evolution 

In this section we define a dynamical model for the temporal evolution of the firms' market 
shares in multiple interacting markets. The model can be regarded as a random element 
whose realizations are right-continuous functions from [0,oo) to the space (X Dn , 3fr Dn ), 
D n G N being (10), and we refer to it as a particle system, since it explicitly models the 
evolution of the share units, or particles, in several markets. For ease of presentation, 
we approach the construction by first considering a single market for a fixed number n 
of share units, and then extend it to a collection of markets. The investigation of the 
asymptotic properties as n — > oo is, instead, the object of Section 4. 

For any fixed n > 1, consider a vector x = x^™) = (x\ , . . . , x n ) £ X™, and let 
(x* , . . . , x* K ) denote the K n < n distinct values in x, with Xj having multiplicity rij. 
The elements of (x^ , . . . , x* K ) represent the K n firms operating in the market at a given 
time. Here x* is a random label to be seen as a unique firm identifier. The vector x 
represents the current market configuration, carrying implicitly the information on the 
shares. Namely, the fraction of elements in x equal to x* is the market share possessed 
by firm j. Here n represents the level of share fractionalization in the market. Dividing 
the market into n fractions is not restrictive, since any share can be approximated by 
means of a sufficiently large n. See Remark 5.1 below for a discussion of the implications 
of this assumption on the computational costs. 

Define now a Markov chain taking values in X™ as follows. At each step an index i 
is chosen from {1, . . . with probability 7^ > for i = 1, . . . ,n, with X)™=i 7»M = ^ - 
Equivalcntly, let 7y(n n ) be the probability that firm x* loses an nth fraction of its 
market share at a certain transition, where 7j(n„) depends on the frequencies n„ = 
(n\, . . . ,riK n )- That is, firm x* undergoes a shock whose probability is idiosyncratic, 
depending on the firm itself and on the current market configuration, summarized by 
the vector of frequencies. Different choices of 7j(n„) reflect different market regulations, 
possibly imposed by the policy maker. We provide some examples: 

(1) 7j(n n ) = 1/ K n : neutrality. All firms have equal probability of undergoing a shock; 

(2) 7j (n„) = rij/n: firms with higher shares are the weakest, with a flattening effect on 
the share distribution. This parametrization is also useful in population genetics 
contexts, where particles represent individuals; 
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(3) 7j(n n ) = (1 — rij/n)/(K n — 1) when K n > 2: firms with higher shares are the 
strongest. The probability of losing shares is decreasing in the firms' positions 
in the market; 

(4) 7j(n n ) = l(maxirii < nC)7 J (n„) + l(rij > nC) for some constant < C < 1, where 
1(A) is the indicator function of the event A. The probability of selecting x* is 
7j(n n ) provided no firm controls more than C% of the market. If firm x* controls 
more than C% of the market, at the following step, x* is selected with probability 
one. Thus C is an upper bound imposed by the policy maker to avoid dominant 
positions. Incidentally, there is a subtler aspect of this mechanism which is worth 
commenting upon. It will be seen later that there is positive probability that 
the same firm acquires the vacant share again, but this only results in picking 
again x* with probability one, until the threshold C is crossed downwards. This 
seemingly anomalous effect can be thought of as the viscosity with which a firm 
in a dominant position gets back to a legitimate status when condemned by the 
antitrust authority, which in no real world occurs instantaneously 

Suppose now Xi = x* has been chosen in x. Once firm x* looses a fraction of its share, 
the next state of the chain is obtained by sampling a new value for Xi from (4), leaving all 
other components unchanged. Hence the zth fraction of share is reallocated, according to 
the predictive distribution of Xi\x(_j\ , either to an existing firm or to a new one entering 
the market. 

Remark 3. 1 . The above Markov chain can also be thought of as generated by a Gibbs 
sampler on q n (dx\ dx n ). This consists of sequentially updating one randomly selected 
component at a time in (x\, . . . ,x n ) according to the component-specific full conditional 
distribution g n .i(dxi|x(_j)). The Gibbs sampler is a special case of a Metropolis-Hastings 
Markov chain Monte Carlo algorithm, and, under some assumptions satisfied within the 
above framework, yields a chain which is reversible with respect to q n (dxi, . . . ,dx ra ), 
hence also stationary. See [18] for details and Appendix A for a brief account. 

Consider now an arbitrary collection of markets, indexed by r, r',r", . . . £ I, so that 
the total size of the system is (10), and extend the construction as follows. At each 
transition, a market r is selected at random with probability g r , and a component of 
(x\, . . . ,x r n ) is selected at random with probability 7^. The next state is obtained by 
setting all components of the system, different from x\, equal to their previous state, and 
by sampling a new value for x\ from (14). Choose now 




(15) 



r'Gl 



where 6 > 0, ir 6 [0, 1], v^, is a non-atomic probability measure on X, 



n 




(16) 
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and m(r,r'):lxl->[0,l] is such that 



m(r, r) = 0, 



\J m(r, r') = 1. 

r'GX 



(17) 



In this case (14) becomes 



(18) 



7 / 



+ 0(1 - 7r)/3„«) ^ m(r,r')^(d<) + ^/3„(4)^(d<) 




with normalizing constant q Dn i = 0(?i) when /?„ = 1 + 0(n 1 ). By inspection of (18), 
there are three possible destinations for the allocation of the vacant share: 

(i) A new firm is created and enters the market. The new value of the location x\ is 
sampled from vq, which is non-atomic, so that x\ has (almost surely) never been 
observed. Here v$ is common to all markets. The possibility of choosing different 
i/o,n t el, is discussed in Section 5 below. 

(ii) A firm operating in the same market r expands its share. The location is sampled 
from the last term, which is a weighted empirical measure of the share distribution 
in market r, obtained by ignoring the vacant share unit x\. 

(iii) A firm operating in another market r' either enters market r or expands its current 
position in r. The location is sampled from the second term. In this case, an index 
r' 7^ r is chosen according to the weights m(r, •); then within r' a firm x T j» is chosen 
according to the weighted empirical measure 



If the cluster associated to x\ has null frequency in the current state, we have an 
entrance from r'; otherwise, we have a consolidation in r of a firm that operates, 
at least, on both those markets. 

We can now provide interpretation for the model parameters: 

(a) 9 governs barriers to entry: the lower the 9, the higher the barriers to entry, both 
for entrance of new firms and for those operating in other markets. 

(b) 7r regulates sunk costs: given 9, a low tt makes expansions from other sectors more 
likely than start-up of new firms, and vice versa. 

(c) m(r, •) allows us to set the costs of expanding to different sectors. For example, 
it might represent costs of technology conversion a firm needs to sustain or some 
regulation constraining its ability to operate in a certain market. Tuning m(-, •) on 
the base of some notion of distance between markets allows us to model these costs, 
so that a low m(r, r') implies, say, that r and r' require very different technologies, 
and vice versa. 



n 



k=l 
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(d) j3 n is probably the most flexible parameter of the model, which, due to the minimal 
assumptions on its functional form (see Section 2), can reflect different features of 
the market, implying several possible interpretations. For example, it might repre- 
sent competitive advantage. Since /3„ assigns different weights to different locations 
of X, the higher /3 n (x*j, ), the more favored is x*j, when competing with the other 
firms in the same market. Here, and later, x T j» denotes the jth firm in market r. It 
is, however, to be noted that setting /3 = 1 does not imply competitive neutrality 
among firms, as the empirical measure implicitly favors those with higher shares. 
More generally, observe that the model allows us to consider a weight function of 
type P n (x r kl ii r ), where \i r is the empirical measure of market r, making /?„ depend 
on the whole current market configuration and on x r k explicitly. This indeed allows 
for multiple interpretations and to arbitrarily set how firms relate to one another 
when competing in the same market. For example, this more general paramctriza- 
tion allows us to model neutrality among firms by setting p n (x k , /j, r ) = with 

being the number of share units possessed by firm j in market r. 

(e) Weights 7 nj i can model barriers to exit, if appropriately tuned (see also points (1) 
to (4) above). For example, setting 7j(n„) very low (null) whenever nj, or nj/n, is 
lower than a given threshold makes the exit of firm xj» very unlikely (impossible) . 

The function /? ra , in point (d) above, will represent the crucial quantity which will 
be used for introducing explicitly the micro-foundation of the model. However, we do 
not pursue this here since we focus on generality and adaptability of the model. The 
micro-foundation will be the object of a subsequent work. 

4. Infinite dimensional properties 

From a qualitative point of view the outlined discrete-time construction would be enough 
for many applications. Indeed Section 5 below presents two algorithms which generate 
realizations of the system and are used to perform a simulation study, based on the 
above description. It is, however, convenient to embed the chain in continuous time, 
which makes the investigation of its properties somewhat simpler and leads to a result 
of independent mathematical interest. This will enable us to show that an appropriate 
transformation of the continuous time chain converges in distribution to a well-known 
class of processes which possess nice sample path properties. To this end, superimpose 
the chain to a Poisson point process with intensity A„, which governs the waiting times 
between points of discontinuity. The following proposition identifies the generator of the 
resulting process under some specific assumptions which will be useful later. Recall that 
the infinitesimal generator of a stochastic process {Z(t),t > 0} on a Banach space L is 
the linear operator A defined by 

Af = \vca- t nf{Z{t))\Z{Q)}-f{Z{Q))] 

with domain given by the subspace of all / £ L, for which the limit exists. In particular, 
the infinitesimal generator carries all the essential information about the process, since 



Bayesian modeling of market dynamics 1 1 

it determines the finite-dimensional distributions. Before stating the result, we need to 
introduce some notation. Let B(K) be the space of bounded measurable functions on X, 
and (g r )rex be a sequence with values in the corresponding simplex 



(19) 



A # z = I (g r )rex- p r > 0, Vr G I, ^ g r = 1 

rex 

Furthermore, let q . be as in (18), with q D ^ . its normalizing constant, and let 

p n (z) = l + tr(z)/n, ctgB(X), (20) 

Cn,r,i = KQr1n,i/q Dn<i - (21) 

Define also the operators 

%(x|z) = (x 1 ,...,x i -i,z,Xi+i,...,x n ), (22) 

M n g(w) = [[g(y)-g(w)}(l + a(y)/n)Mdy), geB(X), (23) 



G n > r g(w)= l\g(y)-g(w)}(l + a(y)/n)n r ,(dy), geB(X), (24) 
and denote by 

r) ri , M?J, G^f'f (25) 

such operators as applied to the ith coordinate of those in x which belong to r. For 
instance, if y = (y\ , 2/2 7 i/3 ? 2/^ )j where 2/2,2/3 belong to market r and the others to r', 
then r)r 2 (y\z)=r) 3 (y\z) = (y[ ,yl,z,yl). 

Proposition 4.1. Let X^ Dn \-) = {X^ Dn '(t),t > 0} be the right- continuous process with 
values in X 15 " which updates one component according to (18) at each point of a Pois- 
son point process with intensity A„. Then X^ Dn \-) has infinitesimal generator, for 
f € B(X D »), given by 



reX K. i=l 



i=l 

n 

9(1 -tt)J2 m(r, r') J2 C n,r,iG n r f /(x) 



i=l 



£ ^[/(^(X|4))-/(X)] 

l<k=£i<n 

- £ a w( r(4)[/(^(xi4))-/(x)]i 

n * — ' 

Kk^i<n ) 



(26) 
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With respect to the market dynamics, generator (26) can be interpreted as follows. The 
first term governs the creation of new firms, obtained by means of operator (23) which 
updates with new values from v$ . The second regulates the entrance of firms from other 
markets, via operator (24) and according to the "distance" kernel m(-,-). The last two 
terms deal with the expansion of firms in the same market. These parallel, respectively, 
points (i), (iii) and (ii) above. 

Consider now the probability-measure- valued system associated with X^ Dn \-), that 
is, F<»(.) = {yW(t),f >0}, where 

r(")(t) = (MO, /v (*),•••) (27) 

and /j, r is as in (16). Y^ n \t) is thus the collection of the empirical measures associated to 
each market, which provides aggregate information on the share distributions at time t. 
The following result identifies the generator of Y^ n \-), for which we need some additional 
notation. Let 

7i[ fc ]=n(ra-l)---(ra-fe+l), n [Q] = l. (28) 

For every sequence (r%, . . . , r m ) £ X m , m 6 N, and given r £ X, define k r = Y^jLi l( r j = r ) 

to be the number of elements in (r±, . . . , r m ) equal to r. Define also pt!^ r ^ and p^ m ^ to be 
the probability measures 



— . )> ( 29 ) 



l<»7-,ir^-"#»T-,fc r <n 



(so) 



and let 



<M,i) = y /d M w, feB(x m ). (3i) 

Finally, denote o rk {-) = o~(x r k ) and 

rk,i = {rk,n), (32) 
with n as in (25), and define the map $ fc?; : B(K n ) ->■ ^(X"- 1 ) by 

$kif{x\, ■ ■ ■ , X n ) = f(xi, . . . , Xi-l, X k , X i+ i, . . . , x n ). (33) 

Proposition 4.2. Let Y {n \-) be as in (27). Then, for <f> m ((j,) as in (31), m < D n and 
under the hypothesis and notation of Proposition 4-1, the generator ofY^( ) is 
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m „ 
r' i=\ 

E C n,r,iJ(^r k J-f)d^ (34) 

"EE^ / ^(-)(*r k , t /-/)dM (m) 
i=l fc^i ^ 



The interpretation of (34) is similar to that of (26), except that (34) operates on the 
product space ^(X)* 1 instead of the product space of particles. Let £? n (X) C £?(X) 
be the set of purely atomic probability measures on X with atom masses proportional 
to n _1 , Z)^)(x)#i([0, oo)) be the space of right-continuous functions with left limits from 
[0, oo) to ^(X)* 1 and C&>(x)*z{[0, oo)) the corresponding subset of continuous functions. 
The following theorem, which is the main result of the section, shows that the measure- 
valued system of Proposition 4.2 converges in distribution to a collection of interacting 
Fleming- Viot processes. These generalize the celebrated class of Fleming- Viot diffusions, 
which take values in the space of probability measures, to a system of dependent diffusion 
processes. See Appendix A for a brief review of the essential features. Here convergence 
in distribution means weak convergence of the sequence of distributions induced for each 
n by F"(-) (as in Proposition 4.2) onto the space D^( X )#i([0, oo)), to that induced 
on the same space by a system of interacting Fleming- Viot diffusions, with the limiting 
measure concentrated on C^a(x)#x([0, oo)). 

Theorem 4.3. Let (•) = {yW (t), t > 0} be as in Proposition with initial distri- 
bution Q n E (^"(X))^ 1 , and let Y(-) = {Y(t),t > 0} be a system of interacting Fleming- 
Viot processes with initial distribution Q € (^(X))^ 1 and generator defined in Ap- 
pendix A by (S8)-(39). Assume X = [0,1], a(-,-) =m(-,-) and M*(x,dy) = u (dy). If 
additionally a in (39) is univariate, X n = 0(n 2 #I) and Q n =>• Q, then 

F(")(.)=>y(.) asn^oo 
in the sense of convergence in distribution in C^(x)#z([0, oo)). 



5. Algorithms and simulation study 

In this section we device suitable simulation schemes for the above constructed systems 
by means of Markov chain Monte Carlo techniques. This allows us to explore different 
economic scenarios and perform sensitivity analysis on the effects of the model parameters 
on the regime changes. Remark 3.1 points out that the discrete representation for a single 
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market can be obtained by means of Gibbs sampling the joint distribution q n ^ in (3). 
A similar statement holds for the particle system in a multi market framework. The 
particle system in Section 3 is such that after a market r and an item x\ are chosen with 
probability g r and 7^ i respectively, a new value for x[ is sampled from 

n 

9 D „,.(d<|2:(-0) oc/3„«)a I( _ r) (da£) + ^TPn(xl)S x r (da£), 

which selects the next ownership of the vacant share, and all other items are left un- 
changed. It is clear that q D i is the full conditional distribution of x\ given the current 
state of the system. Since the markets, and the particles within the markets, are updated 
in random order, it follows immediately that the particle system is reversible, hence 
stationary, with respect to (11). 

Algorithm 1 is the random scan Gibbs sampler which generates a sample path of the 
particle system with the desired number of markets. Here we restrict to the case of a = 0, 
which implies that the normalizing constant q D i is 6 + n — 1. 



Algorithm 1 

Initialize; then: 

1. select a market r with probability g r ; 

2. within r, select x\ with probability 7^ 

3. sample u ~ Unif(0, 1); 

4. update x\: 

a. if u < e _^_ 1 , sample x\ ~ v a ; 

b. if u > gr~3 ) sample uniformly an x r k , k^i, within market r and set x\ = x r k ; 

c. else: 

i. select a market r' with probability m(r, r'); 

ii. sample uniformly an x'j within market r' and set x\ = x r j ; 

5. go back to 1. 



Remark 5. 1 . Note that the fact that updating the whole vector implies sampling from 
n different distributions does not lead to an increase in computational costs if one wants 
to simulate from the model. Indeed, acceleration methods such as those illustrated in [25] 
can be easily applied to the present framework. 

As previously mentioned, setting a = 0, hence /3 = 1, as in Algorithm 1, does not lead 
to neutrality among firms, determining instead a competitive advantage of the largest 
(in terms of shares) on the smallest. A different choice for /3 allows us to correct or 
change arbitrarily this feature. For example, choosing /3(a^„ ,fi r ) — nj 1 , where rij is the 
absolute frequency associated with cluster Xj, , yields actual neutrality. Observe also that 
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sampling from (18), which is composed of three additive terms, is equivalent to sampling 
either from 



with probability 



from 



with probability 



or from 



with probability 



/3n(2/,MrH(dy), 



Er'Gl TO ( r ' r E"=l Aifaj' , /Xr) 

Efc^^(4^r)^(dO 
1 ™ 



(35) 



(36) 



(37) 



where the normalizing constant g Dn . is given by 

8tt [ p n (x)v (dx) + e(l-ir)^2m(ry) f /3„(a:)/v (dx) + ^ /U4)- 

Once the functional forms for f3 and m are chosen, computing q Dn . is quite straightfor- 
ward. If, for example, X = [0, 1], and the type of an individual admits also interpretation 
as index of relative advantage, then one can set (3(x) = x, and q D . becomes 

n 

0nD o + 0(1 - tt) m (r, r')x r ' + ^ x r k , 

r'el k^i 

where is the mean of and x r is the average of the components of market r' . To this 
end, note also that the assumption of Vq being non-atomic can be relaxed simplifying the 
computation. Algorithm 2 is the extended algorithm for f3 n ^ 1 . 

In the following we illustrate how the above algorithms produce different scenarios 
where economic regime transitions are caused or affected by the choice of parameters, 
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Algorithm 2 

Initialize; then: 

1. select a market r with probability g r ; 

2. within r, select x\ with probability 7^ 

3. sample u ~ Unif (0, 1); 

4. update x\: 

a. if u < q^ 1 .irO J f3 n {y, A t r)' y o(dj/), sample x\ from (35); 

b. if u > 1 - Sfc 7 ti^n(^fe)A 4 '-)) sample^ from (37); 

c. else sample a;[ from (36); 

5. go back to 1. 



which can be structural or imposed by the policy maker during the observation period. 
We first consider a single market and then two interacting markets, and for simplicity 
we confine to the use of Algorithm 1. As a common setting to all examples we take 
X= [0,1], n = 500, Uq to be the probability distribution corresponding to a Beta(a,6) 
random variable, with a,b > 0, with the state space discretized into 15 equally spaced 
intervals. The number of iterations is 5 x 10 5 , of which about 150 are retained at increasing 
distance. Every figure below shows the time evolution of the empirical measure of the 
market, which describes the concentration of market shares, where time is in log scale. 

Figure 1 shows a single market which is in an initial state of balanced competition 
among firms, which have similar sizes and market shares: this can be seen by the flat 




Figure 1. High sunk costs progressively transform a perfectly competitive market into an 
oligopoly and then into a monopoly. 
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Figure 2. An Oligopoly becomes a competitive market after the policy maker reforms the sector 
regulation (threshold 1), and concentrates again after the reform is abolished (threshold 2). 

side closest to the reader. As time passes, though, the high level of sunk costs, deter- 
mined by setting a low 8, is such that exits from the market are not compensated by 
the entrance of new firms, and a progressive concentration occurs. The competitive mar- 
ket first becomes an oligopoly, shared by no more than three or four competitors, and 
eventually a monopoly. Here uq corresponds to a Beta(l,l) and 6 = 1. The fact that 
the figure shows the market attaining monopoly and staying there for a time greater 
than zero could be interpreted as conflicting with the diffusive nature of the process 
with positive (although small) entrance rate of new firms (mutation rate in population 
genetics terms). In this respect it is to be kept in mind, as already mentioned, that 
the figure is based on observations farther and farther apart in time. So the picture 
does not rule out the possibility of having small temporary deviations from the seem- 
ing fixation at monopoly, which, however, do not alter the long-run overall qualitative 
behavior. 

In Figure 2 we observe a different type of transition. We initially have an oligopolis- 
tic market with three actors. The structural features of the market are such that the 
configuration is initially stable, until the policy maker, in correspondence to the first 
black solid line, introduces some new regulation which abates sunk costs or barriers to 
entry. Note that in the single market case the parameter 8 can represent both, since 
this corresponds to setting n = 1 in (15), while in a multiple market framework we can 
distinguish the two effects by means of the joint use of 8 and tt. Here all parameters are 
as in Figure 1, except 8, which is set equal to 1 up to iteration 200, equal to 100 up to 
iteration 4.5 x 10 4 and then equal to 0. The concentration level progressively decreases 
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and the oligopoly becomes a competitive market with multiple actors. In correspondence 
of the second threshold, namely the second black solid line, there is a second regulation 
change in the opposite direction. The market concentrates again, and, from this point 
onward, we observe a dynamic similar to Figure 1 (recall that time is in log scale, so 
graphics are compressed toward the farthest side). The two thresholds can represent, for 
example, the effects of government alternation when opposite parties have very different 
political views about a certain sector. 

We now proceed to illustrate some effects of the interaction between two markets with 
different structural properties and regulations when some of these parameters change. 
Figure 3 shows three scenarios regarding a monopolistic (left) and a competitive market 
(right). In all three cases v§ corresponds to a Beta(l, 1) for both markets. Case 1 repre- 
sents independent markets, due to very high technological conversion costs or barriers to 
entry, which is for comparison purposes. Here 8 a =Q,6b = 100 and 7Tb = 1. In Case 2 the 
monopolistic market has low barriers to entry, while (2b) is still closed, and a transition 
from monopoly to competition occurs. Here 8 a = 30,6^ = 100, 7r a = 0.01, lib = 1- Case 3 
shows the opposite setting, that is, a natural monopoly and a competitive market with 
low barriers to entry. The monopolist enters market (3b) and quickly assumes a dominant 
position. Here 9 a = : b = 100, 7Tb = 0.7. Recall, in this respect, the implicit effect due to 
setting P = 1, commented upon above. 

Case (2a) in Figure 3 suggests another point. The construction of the particle system 
by means of the hierarchical models defined in Section 2 compels us to have the same 
centering measure i>o, which generates new firms for all markets. In particular, this makes 
it essentially impossible to establish, by mere inspection of Figure 3 (2a), whether the 
transition is due to new firms or to entrances from (2b) . Relaxing this assumption on vq 
partially invalidates the underlying framework above, in particular, due to the fact that 
one loses the symmetry implied by Lemma 2.1. Nonetheless the validity of the particle 
system is untouched, in that the conditional distributions of type (14) are still available, 
where now v$ >r , i n place of a common vq, is indexed by r G X. This enables us to appreciate 
the difference between the two above mentioned effects. If one is willing to give a specific 
meaning to the location of the point ieX which labels the firm, then v$ :r ^ i^.r' can 
model the fact that, say, in two different sectors, firms arc polarized on opposite sides of X, 
which, in turn, represents some measurement of a certain exogenous feature possessed by 
those firms. Consider a monopoly and a competitive market, where we now take i^o, a and 
^o,f, to be the probability measures corresponding to a Beta(2, 4) and a Beta(4, 2) random 
variable for the monopolistic and competitive market, respectively. We are assuming that 
firms on the left half of the state space have a certain degree of difference, with respect to 
those on the other side, in terms of a certain characteristic. Figure 4 shows the different 
impact of barriers to entry and sunk costs on the monopolistic market, due to the joint 
use of 7r and 6, thus splitting Figure 3(2a) into two different scenarios. The competitive 
market is composed by firms which are polarized toward the right half of the state space, 
meaning, for example, that they have a high level of a certain feature. Then case 1 of 
Figure 4 shows the monopoly when sunk costs are high, but barriers to entry are low, 
so that the concentration is lowered by entrance of firms from the other market, rather 
than from the creation of new firms from within; while case 2 shows the effects of high 
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Figure 3. Effects of parameters' change in interacting monopolistic and competitive markets, 
(la) and (lb) are both closed, hence independent, markets. (2b) is closed, but (2a) has low 
barriers to entry (n « 0), and firms from (2b) progressively lower the concentration in (2a). 
(3b) has low barriers to entry, so that the monopolist of (3a) enters the market and conquers a 
dominant position. 

barriers to entry and low sunk costs, so that a transition to a competitive regime occurs 
independently of (2b). The parameters for case 1 are 9 a = 30, 9b = 100, 7r a = 0, tvi, = 1, 
while for case 2 we have 8 a = 30, 9b = 100, 7r a = 1, 7T& = 1. 
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Figure 4. Firms in the competitive market (right) are polarized towards the right half of the 
state space, (la) is a monopoly with high sunk costs and low barriers to entry, so firms from 
(lb) enters market (la). (2a) is a monopoly with high barriers to entry and low sunk costs, so 
that a transition to a competitive regime occurs independently of (2b). 

6. Concluding remarks 

In this paper we propose a model for market share dynamics which is both well founded, 
from a theoretical point of view, and easy to implement, from a practical point of view. In 
illustrating its features we focus on the impact of changes in market characteristics on the 
behaviors of individual firms taking a macroeconomic perspective. An enrichment of the 
model could be achieved by incorporating exogenous information via sets of covariatcs. 
This can be done, for example, by suitably adapting the approach recently undertaken 
in [34] to the present framework. Alternatively, and from an economic viewpoint, more 
interestingly, one could modify the model adding a microeconomic understructure: this 
would consist of modeling explicitly the individual behavior by appropriately specifying 
the function /3 n at Point (d) in Section 3, which can account for any desired behav- 
ioral pattern of a single firm depending endogenously on both the status of all other 
firms and the market characteristics. This additional layer would provide a completely 
explicit micro-foundation of the model, allowing us to study the effect of richer types 
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of heterogeneous individual decisions on industry and macroeconomic dynamics through 
comparative statics and dynamic sensitivity analysis. These issues of more economic 
flavor will be the focus of a forthcoming work. 

Appendix A: Background material 

Basic elements on the Gibbs sampler 

The Gibbs sampler is a special case of the Metropolis-Hastings algorithm, which, in turn, 
belongs to the class of Markov chain Monte Carlo procedures; see, for example, [18]. These 
are often applied to solve integration and optimization problems in large dimensional 
spaces. Suppose the integral of /:X — s-R d with respect to ir <G ^(X) is to be evaluated, 
and Monte Carlo integration turns out to be unfeasible. Markov chain Monte Carlo 
methods provide a way of constructing a stationary Markov chain with ir as the invariant 
measure. One can then run the chain, discard the first, say, N iterations, and regard 
the successive output from the chain as approximate correlated samples from n, which 
are then used to approximate f fdx. The construction of a Gibbs sampler is as follows. 
Consider a law n = 7r(dxi, . . . , dx n ) defined on (X n , J£T n ), and assume that the conditional 
distributions 

7r(dxi |xi , . . . , , j . . . , x n ) 

are available for every 1 < i < n. Then, given an initial set of values (x\, . . . , x° ), update 
iteratively 

x\ ~7r(dx]>2,...,x°), 
x\ ~7r(dx 2 |xi,X3,...,x°), 

x n ~ 7r(dx n |x^, . . . , x n _ 1 ), 
x\ ~7r(dxi|z2,...,a£), 

and so on. Under mild conditions, this routine produces a Markov chain with equilib- 
rium law 7r(dxi, . . . ,dx„). The above updating rule is known as a deterministic scan. If 
instead the components arc updated in a random order, called random scan, one also 
gets reversibility with respect to tt. 

Basic elements on Fleming— Viot processes 

Fleming- Viot processes, introduced in [17], constitute, together with Dawson- Watanabe 
superprocesses, one of the two most studied classes of probability-measure-valued diffu- 
sions, that is, diffusion processes which take values on the space of probability measures. 
A review can be found in [15]. 
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A Fleming- Viot process can be seen as a generalization of the neutral diffusion model. 
This describes the evolution of a vector z = {z%)i^s representing the relative frequencies 
of individual types in an infinite population, where each type is identified by a point in 
a space S. The process takes values on the simplex 



[ (4 eS 6 [0,l] s : z i >0,J2zi = ~ i ) 



and is characterized by the infinitesimal operator 



5 2>fo-*,) g^ + E^, 



defined, for example, on the set C(S) of continuous functions on S, if S is compact. Here 
the first term drives the random genetic drift, which is the diffusive part of the process, 
and bi(z) determines the drift component, with 

where <?y is the intensity of a mutation from type i to type j and cr^ = o~ji is the selection 
term in a diploid model. This specification is valid for S finite, which yields the classi- 
cal Wright -Fisher diffusion, or countably infinite; sec, for example, [13]. Fleming and 
Viot [17] generalized to the case of an uncountable type space S by characterizing the 
corresponding process, which takes values in the space &(S) of Borel probability mea- 
sures on S, endowed with the topology of weak convergence. Its generator on functions 
<Mm) = F((fi,fJ,),. . . , (f m , m)) = F((f , /it)), where F E C 2 (R m ), fx,..., f m continuous on 
S and vanishing at infinity, for m > 1, and (/, fi) = f fd/i, can be written 

i,J=l 

m rn 

+ ^ ( A//, , fi) F Zi ( (f , M ) ) + ]T ( ( (/, o 7r)a, /i 2 ) - </« , /i) (a, M 2 ) ) F Zi ( (f , /i> ) , 

i=l i=l 

where fj, 2 denotes product measure, ir is the projection onto the first coordinate, M is the 
generator of a Markov process on S, known as the mutation operator, a is a non-negative, 
bounded, symmetric, Borel measurable functions on S 2 , called selection intensity function 
and F Zi is the derivative of F with respect to its ith argument. Recombination can also 
be included in the model. 



Interacting Fleming— Viot processes 



Introduced by [40], and further investigated by [7] and [6], a system of interacting 
Fleming-Viot processes extends a Fleming-Viot process to a collection of dependent 
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diffusions of Fleming- Viot type, whose interaction is modeled as migration of individu- 
als between subdivided populations. Following [6], the model without recombination can 
be described as follows. Let the type space be the interval [0, 1]. Each component of the 
system is an element of the set &([0, 1]), denoted fi r and indexed by a countable set X 
of elements r, r' , . . . . For F : (^([0, l])) x -> K of the form 

F (l J -)= ■■ f(xi,... 1 x m )fi ri (dxi)---fj Tm (dx m ) (38) 

J [0,1] J [0,1] 

with / 6 C([0, l] m ), (ri, . . . ,r m ) S (X) m , m G N, the generator of a countable system of 
interacting Fleming- Viot processes is 



9F{n), , , dF(u), ' 

— ^ y )M* (x,dy) - -^(x) 



H r (dx) 



c V a(r,r') / (/x r 



2 



Mr )(dx)^M(x) 

OfJL r 



^ {x,y)Q llr {Ax,dy) 



[0,1] J [0,1] d^ r dn r 

dF(fx) 



(39) 



(z)o-(y, z)ii r (dy)Q llr {dx, dz) 

[0,1] J [0,1] J [0,1] Cf^r 

where the term Q flr (dx,dy) = fj, r (dx)S x (dy) — fj, r (dx)fj, r (dy) drives genetic drift, M*(x,dy) 
is a transition density on [0, 1] x &([0, 1]) modeling mutation, &([0, 1]) is the Borel sigma 
algebra on [0, 1], a(-, •) on X x X such that a(r, r') <E [0, 1] and ^ r a(r, r 1 ) — 1 is a transition 
kernel modeling migration and cr(-, •) is a bounded symmetric selection intensity function 
on [0, l] 2 . The non-negative reals q,c,d,s represent, respectively, the rate of mutation, 
immigration, resampling and selection. Let the mutation operator be 



Mf(z) = J if(y) - f(z)]M*(x, dy), fe B{X) (40) 
and the migration operator be 

G r 'f(z)= [{f(y)-f(z)]fi r ,(dy), /eB(X), (41) 
for r' <E X. Using this notation, and when F is as in (38), (39) can be written 

GF(p) = I <lf2 l ■ f M rjdfl ri ■ ■ ■ d M r m 

ren N { i= i-W] J [OA] 

+ c V a(r,r')V / ••• / / dfi ri ■ ■ ■ dfi rm 

r'en N 1=1 Ml] M,i] 
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mm* „ 

+ d EE / ■'•/ ^r k J-f)d^ ri ---d^ rm (42) 

t =lk^J[OA] J [OA] 

m „ ,. 

+ S E / " / (°n,m+i (■,■)/ 

4=1 J[0,1] J[0,1] 

-t7 m+ l )TO+2 (-,-)/)d/X ri ••• dn rm dn r dfX r 1 , 

where Mj and GJ are M and G r applied to the jth coordinate of /, r, is as in Propo- 
sition 4.1, rfc.i as in (32) and as in (33). When I is single- valued, (42) simplifies 
to 

m mm 

GFQt) = qJ2(Mif,» m ) + - /./O 

i=l i— 1 fcy^i 

m 

+ S ^((a i>m+1 (-, -)/,^ i+1 ) - (*m+l,Tn+2(-, 0/^ m+2 ». 

which is the generator of a Fleming- Viot process with selection with F(/j,) = (f,/j, m ), 

/eC([0,l]). 

Appendix B: Proofs 



Proof of Proposition 4.1. The infinitesimal generator of the X™-valued process de- 
scribed at the beginning of Section 3 can be written, for any / £ £?(X n ), as 

n ~ 

4/(x) = X n ^2ln,i / [/(%(x|y)) -/(x)]g„ i i(d?/|x ( _ i )), (43) 

z=l 

where g rii i(dj/|x(_ i )) is (4) and rji is as in (22). Within the multi-market framework, (43) 
is the generator of the process for the configuration of market r, say, conditionally on all 
markets r' £ X, r' ^ r, and can be written 

A D J(x r \I(-r)) = \ n J2Yn,i] [f(Vi(x r \y))-f(x r )]qn n Ady\l(- X r)), (44) 

i—l 

where I(—r) and I(-x^) are as in (12) and (13), r f ni are the market-specific removal 
probabilities and q Dn 4 (dj/|X(— x\)) is (14). Then the generator for the whole particle 
system, for every / £ B(X Dn ), is 

A D J(^ = X n ^Qr^n,i / [/(^(x| 2 /))-/(x)]^ i (dy|X(-<)), (45) 
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where r\ ri is as in (25). Setting now f3 n as in (20), (45) becomes 



rei I i=i \ n 



E a i ,,,,[/(^(x|4))-/(x)] (46) 

l<fc^i<ri 

- E C„, r>i( 7K)[/(r ?n (x|^))-/(x)]} J 
n * — ' 

l<k=£i<n ) 



with C„ >ri i as in (21). Substituting (15) in (46) yields 

A D J(x) = eI^E^ / [/(^(x|y)) - /(x)] (l + ^f) 



^o(dy) 



0(1 - tt) E »n(r, r') E C »rt J Iffa* ^ ~ ^ 



x N+-M ) AV (dy) 



+ E C„ )r>i [/(77 ri (x|4))-/(x)] 

l<k^i<n 

+ \ E a w a(xD[/(^(x|4))-/(x)]l. 

By means of (23) and (24), with Mf / and G"' r '/ denoting, respectively, M" and G"' r ', 
applied to the ith coordinate of /, and M™/ and G™.' r / interpreted according to (25), 
(47) can be written 



rei (. 1=1 



A D J(X) = E<| ^ECn,r,iM r y(x) 
i=l 

n 

0(1 - tt) E m(r, r') E ^r,^' 7(x) 



i=i 



+ E ^r,i[/(»7n(x|4))-/(x)] 

l<fc^i<n 

+ - E c-„,^(4)[/(^(x|4))-/(x)]l. 

Kk^i<n ) 



□ 
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Proof of Proposition 4.2. For k <n, let nm be as in (28), and define the probability 
measure 



/' 



CDh) _ 



reX 11 l<i r .i#---#ir,fc<n 



^( x i r . 1 --- a: i r fc )' 



(48) 



where Dk is as in (10). Define also 

4> D M = (f,» {Dk) ), feB{x D ») 

and 

A D ^ Dk (^) = (A D J,^) 1 



(49) 



where (/, — J f dfi. Then Ajj^cpn^ (n) is the generator of the (^(X)) ^-valued system 
(27), which from (26), letting / e B(X D ") in (49), can be written 



A Dn cj) Dn ((i) = £ 



rex 



(50) 



8irJ2C n , r ,(M? l f^ {Dn) ) 

n 

+ 9(l-ir)Y / Hry)J2 C ^AGr; r 'f^ {D ' l) ) 

r' i—1 

+ £ C n ,rA®r k J - f,V iDn) ) 
\<k^i<n 

n £ — ' 

l<fe/i<n 

where a rk (■) denotes c( x fc) and <£>fci is as in (33). Note now that for / e £?(X TO ), m < £>„ 
we have 



M"/ = /, G£ r / = /, $r k J = f, ifi>m 



and 



Given (29) and (30), it follows that when / G ^(X" 1 ), m < L>„, (50) can be written 

f rn 

A Dn </> m (lJ.) = £<j ^J2 Cn ^ M rJ^ {m) ) 
i=l 

m 

9(1 - tt) £ a (r, r') £ C^G^f, M (m) > 



(=1 



£ Gn.rA^f-f,^) 
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m k r 



n 



X] Cn,r,i(o- m+ l(-)(*TO+l,n/ - /),M (m) Mr) 



□ 



i=l J 

Proof of Theorem 4.3. For / e B(X k ), k>l, let ||/|| = su P;! , eXfc \f(x)\. Observe that 
(23) and (24) converge uniformly, respectively to (40) and (41), as n tends to infinity, 
implying 

\\(M?J,l* lm) ) - (MrJ^ (m) )\\ -> , / G B(X m ), 

||<G™f7,M (m) ) ~ (</,/i (m) )||^0, /eB(r). 

Here the supremum norm is intended with respect to the vector x S X m of atoms in // m ) , 
with as in (30). Let now fif^ be as in (29), so that jj, r = nT 1 Y^i=i Then it is 
easy to check that 

\\{f,ti**))-{f,tfr)\\^o, /eS(X fc -), 

as n — > oo, where /i r denotes a fc r -fold product measure /j, r x ■ • • x fj, r , and that 

IK/,/i (m) >-</,/i xm )||->0, /ei?(x m ), 

as n — > oo , where we have denoted 

We also have, from (21) C n , ri j = \ n Qrln i/^o^ n where A„ is the Poisson rate driving the 
holding times, g r = 0(#I _1 ) and 7^ i = OirT 1 } are the probability of choosing market 
r and x^ respectively during the update, and q D . = 0(n) is the normalizing constant of 
(18). Then choosing A„ = 0(nD n ) = 0(n 2 #Z) implies C„ r j — > 1 as n — > 00. Finally, let 

^ ra GB(#(X m )) be 



(Pm(v) = (/,M xm ) =/ ••• / /(xi, . . . ,x m )// ri (dxi) • • -/x rm (dx m ) (51) 

•'[0,1] ^ [0,1] 

for any sequence (r\, . . . ,r m ) £l m . Then it can be checked that (34) converges, as n 
tends to infinity, to 



.00 = E 



enJ^iMrJ,!**™) +9(1 -tt) E m(r,r')E(</'M Xro ) 



£ ($ rM / - /, + ^(cr ri (•)/ - <r m+ i(-)/, M x "> r ) 

l<k^i<m i—1 



2S 
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which, in turn, implies 

||AD„</> m (ju) - &tp m {ii)\\ — >0 asn^oo. 
Using (51), and letting X = [0, 1], A<p m (/i) can be written 



0ttS2 ■ ■ M r Jdfi ri ■ ■ ■ d^ rm 

~(J[0,l] J[0,1] 

m p p 

+ e(l-ir)y2m(r,r')y2 "• </d^ ••• d^. m 
r'£X l =i- / [o,i] J [0,1] 

(52) 

/ ■•■/ (^ M /-/)cW--d/i rn , 



l<fc^i<m 



Y " / ( £r r i (-)/ _0 m+x(-)/) d A t ri ••• d/i rm d/i r 

.•_ 1 J[0.1l Jfo.il 



i=1 -/[0,l] ^[0,1 

which equals (42) for appropriate values of q, c, d, s and for univariate a. The statement 
with C^>(x)#x([0, oo)) replaced by Dgs^#z([0, oo)) now follows from Theorems 1.6.1 
and 4.2.11 of [14], which, respectively, imply the strong convergence of the corresponding 
semigroups and the weak convergence of the law of Y^ n \-) to that of Y(-). Replacing 
Z)^( X )#t([0, oo)) with C^s( X )#i([0, oo)) follows from [1], Section 18, by relativization of 
the Skorohod topology to C^(x)#x([0, oo)). □ 
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